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We report on some recent results concerning the dynamics of Bose-Einstein condensates, ob- 
tained in a series of joint papers [3 [6] with L. Erdos and H.-T. Yau. Starting from many body 
quantum dynamics, we present a rigorous derivation of a cubic nonhnear Schrodinger equation 
•/^ ' known as the Gross-Pitaevskii equation for the time evolution of the condensate wave function. 



p ,! 1 Introduction 



Bosonic systems at very low temperature are characterized by the fact that a macroscopic fraction 
of the particles collapses into a single one-particle state. Although this phenomenon, known as 
Bose-Einstein condensation, was already predicted in the early days of quantum mechanics, the first 
empirical evidence for its existence was only obtained in 1995, in experiments performed by groups 
led by Cornell and Wieman at the University of Colorado at Boulder and by Ketterle at MIT (see 
ff^ ' [21 HI)- In these important experiments, atomic gases were initially trapped by magnetic fields and 
cooled down at very low temperatures. Then the magnetic traps were switched off and the consequent 
time evolution of the gas was observed; for sufficiently small temperatures, the particles remained 
close together and the gas moved as a single particle, a clear sign for the existence of condensation. 
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, In the last years important progress has also been achieved in the theoretical understanding of 

Q I Bose-Einstein condensation. In [10], Lieb, Yngvason, and Seiringer considered a trapped Bose gas 
consisting of N three-dimensional particles described by the Hamiltonian 



X 



N N 
j=l i<j 

where 14xt is an external confining potential and Va{x) is a repulsive interaction potential with 
scattering length a (here and in the rest of the paper we use the notation Vj = V^j and Aj = Axj). 
Letting N ^ oo and a ^ with Na = ag fixed, they showed that the ground state energy E{N) of 
(|l.ip divided by the number of particle N converges to 

E(N) . ^ , , 

lim — -— = mm tGp(¥') 

N-*oo, Na=ao N <^gL2{r3). ||<^||=i 

where Sqp is the Gross-Pitaevskii energy functional 

^GP(V?) = j^^ (|V(^(x)|2 + Uext(x)|^(x)|2 + 47raol^(x)|^) . (1.2) 
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Later, in [9], Lieb and Seiringer also proved that trapped Bose gases characterized by the Gross- 
Pitaevskii scahng Na = oq = const exhibit Bose-Einstein condensation in the ground state. More 
precisely, they showed that, if V'Af is the ground state wave function of the Hamiltonian (jl.ip and 
if denotes the corresponding one-particle marginal (defined as the partial trace of the density 

matrix = \ipj\[){ipN\ over the last — 1 particles, with the convention that TV 7^'' = 1 for all N), 
then 

^ \4>GP){4>GP\ asiV^oo. (1.3) 

Here cpcp e ^^(M^) is the minimizer of the Gross-Pitaevskii energy functional (jl.2p . The interpre- 
tation of this result is straightforward; in the limit of large N, all particles, apart from a fraction 
vanishing as ^ 00, are in the same one-particle state described by the wave-function cpQp £ L^(M^). 
In this sense the ground state of (jl.ip exhibits complete Bose-Einstein condensation into (/)gp- 

In joint works with L. Erdos and H.-T. Yau (see [51 El E|), we prove that the Gross-Pitaevskii 
theory can also be used to describe the dynamics of Bose-Einstein condensates. In the Gross- 
Pitaevskii scaling (characterized by the fact that the scattering length of the interaction potential is 
of the order 1/N) we show, under some conditions on the interaction potential and on the initial A^- 
particle wave function, that complete Bose-Einstein condensation is preserved by the time evolution. 
Moreover we prove that the dynamics of the condensate wave function is governed by the time- 
dependent Gross-Pitaevskii equation associated with the energy functional ()1.2p . 

As an example, consider the experimental set-up described above, where the dynamics of an 
initially confined gas is observed after removing the traps. Mathematically, the trapped gas can be 
described by the Hamiltonian (jl.ip . where the confining potential I4xt models the magnetic traps. 
When cooled down at very low temperatures, the system essentially relaxes to the ground state ^pN 
of (jl.ip : from [9j it follows that at time t = 0, immediately before switching off the traps, the system 
exhibits complete Bose-Einstein condensation into 0gp in the sense (jl.Sp . At time t = the traps 
are turned off, and one observes the evolution of the system generated by the translation invariant 
Hamiltonian 

N N 
Hn = -^Aj + Va{Xi - Xj) . 
j=l i<j 

Our results (stated in more details in Section [3] below) imply that, if 'ipN,t = e~*^^*VAf is the time 
evolution of the initial wave function V'iv and if 7^ j denotes the one-particle marginal associated 
with ipN,t, then, for any fixed time t € M, 

7/^,1 \y^t){y^t\ as N ^ 00 
where ipt is the solution of the nonlinear time-dependent Gross-Pitaevskii equation 

idtift = -Aipt + 8TTao\ipt\'^ft (1-4) 

with the initial data ^t=o = (pGP- In other words, we prove that at arbitrary time t G M, the 
system still exhibits complete condensation, and the time-evolution of the condensate wave function 
is determined by the Gross-Pitaevskii equation (jl.4p . 

The goal of this manuscript is to illustrate the main ideas of the proof of the results obtained in 
[SlEllT]. The paper is organized as follows. In Section [2] we define the model more precisely, and 
we give a heuristic argument to explain the emergence of the Gross-Pitaevskii equation (jl.4p . In 
Section [3] we present our main results. In Section U] we illustrate the general strategy used to prove 
the main results and, finally, in Sections [5] and [6] we discuss the two most important parts of the 
proof in some more details. 
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2 Heuristic Derivation of the Gross-Pitaevskii Equation 

To describe the interaction among the particles we choose a positive, spherical symmetric, compactly 
supported, smooth function V{x). We denote the scattering length of V by cq. 

Recall that the scattering length of V is defined by the spherical symmetric solution to the zero 
energy equation 

-A + ^V{x)^ f{x) = f{x)-^l as|a;|^cx). (2.1) 
The scattering length of V is defined then by 

ao = lim |x| — \x\f{x) . 



x^oo 



This limit can be proven to exist if V decays sufficiently fast at infinity. Note that, since we assumed 
V to have compact support, we have 

m = 1 - (2-2) 

\x\ 

for |x| sufficiently large. Another equivalent characterization of the scattering length is given by 

8TTao = J dxV{x)f{x). (2.3) 

To recover the Gross-Pitaevskii scaling, we define Viy{x) = N'^V{Nx). By scaling it is clear that 
the scattering length of Vjy equals a = a^/N . In fact if f{x) is the solution to (|2.ip . it is clear that 
fN{x) = f{Nx) solves 

A + Ivn{x)] fN{x)=0 (2.4) 



2 

with the boundary condition /Ar(x) ^ 1 as \x\ oo. From (|2.2p . we obtain 

/at x) = 1 - — — = 1 - -r-r 

iV|x| \x\ 

for \x\ large enough. In particular the scattering length a of is given by a = a^/N . 
We consider the dynamics generated by the translation invariant Hamiltonian 

AT AT 

//7V = ^-A, + ^y7v(xi-x,) (2.5) 

j=l i<j 



acting on the Hilbert space L^(M'^ , dxi . . . dxAr), the bosonic subspace of L (R'^ , dxi . . . dx^r) con- 
sisting of all permutation symmetric functions (although it is possible to extend our analysis to 
include an external potential, to keep the discussion as simple as possible we only consider the 
translation invariant case (j2.5p ). We consider solutions ipN,t of the A^-body Schrodinger equation 

idttpN,t = HNtpN,t ■ (2.6) 

Let ^N,t = \'4^ N ,t) {i' N ,t\ denote the density matrix associated with V'Af.ii defined as the orthogonal 
projection onto ipN,t- In order to study the limit ^ cxd, we introduce the marginal densities of 
^N,t- For A; = 1, . . . , A^, we define the fc-particle density matrix f associated with ipN,t by taking 
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(k) 

the partial trace of 'jN^t over the last N — k particles. In other words, 7]y ^ is defined as the positive 



trace class operator on Lg{Mr ) with kernel given by 

7^1(xfc;x'^)= I dyi]si-ki^N,t{^k,y^N~k)'^N,t{^k^^N-k) ■ (2.7) 



Here and in the rest of the paper we use the notation x = (xi, X2, . . . , xat), = (xi, X2, . . . , x^), 
= (a^'i, x'2, . . . , x'^,), and yiN-k = {xk+i, Xk+2, • • • , xn)- 

We consider initial wave functions ipN,o exhibiting complete condensation in a one-particle state 
ip. Thus at time t = 0, we assume that 

1^nI^\^){v\ asTV^cx). (2.8) 

It turns out that the last equation immediately implies that 

asiV^oo (2.9) 

for every fixed k €N (the argument, due to Lieb and Seiringer, can be found in [9j, after Theorem 1). 
It is also interesting to notice that the convergence (j2.8p (and (12. 9p ) in the trace class norm is 
equivalent to the convergence in the weak* topology defined on the space of trace class operators on 
M or for 1^ ) ; we thank A. Michelangeli for pointing out this fact to us (the proof is based 
on general arguments, such as Griimm's Convergence Theorem). 

Starting from the Schrodinger equation (j2.6p for the wave function ipN,t, we can derive evolution 

(k) 

equations for the marginal densities Jj^ ^- The dynamics of the marginals is governed by a hierarchy 
of coupled equations usually known as the BBGKY hierarchy. 

N k 



j=l i<j 



Xi Xj),^j^, 



(2.10) 



+ {N -k)^ Trfc+i VN{xj - Xfc+i), 7^'^| 
i=i 

Here Tr^+i denotes the partial trace over the {k + l)-th particle. 

(k) 

Next we study the limit A — > 00 of the density 7]y ^ for fixed A; S N. For simplicity we fix /c = 1. 
From (|2.10p , the evolution equation for the one-particle density matrix, written in terms of its kernel 

7^^j(xi;x'i) is given by 



+ (A - 1) / dx2 (VAr(xi - X2) - Vn{x'i - X2)) 7]^^^^ (xi , X2 ; x'^ , X2) . 



(2), , . (2-11) 

N,t' 

Suppose now that 7^^ and 7^\ are limit points (with respect to the weak* topology) of 7^ j and, 
respectively, 7]y ^ as A ^ 00. Since, formally, 

(A - l)F^(x) = (A - l)A2y(Ax) ~ AV(Ax) ^ 60^(2;) with 60 = / dx V{x) 
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(2.13) 



(2.14) 



as A'^ ^ oo, we could naively expect the limit points and to satisfy the limiting equation 

«.,S.(.. - + a;) ,S.(., + / d., - - - ■fS^.^.M ■ 

(2.12) 

From (|2.9p we have, at time i = 0, 

7So(a:i;a;;) = (/j(xi)^(a:;) 

7oo,o(2^i'^2;a;'i,4) = <^(a:^i)v'(a^2)^(x'i)^(4)- 
If condensation is really preserved by the time evolution, also at time t 7^ we have 

Inserting (|2.14p in (j2.12p . we obtain the self-consistent equation 

idt'ft = -A(/?t + 6o|v?ip(/?t (2-15) 

for the condensate wave function ipf This equation has the same form as the time-dependent Gross- 
Pitaevskii equation (jl.4p . but a different coefficient in front of the nonlinearity (60 instead of 87rao). 

The reason why we obtain the wrong coupling constant in (|2.15p is that going from ()2.1ip to 
(|2.12p . we took the two limits 

{N-l)VN{x)^bo6{x) and 7^1 ^ l!^]t (2-16) 

independently from each other. However, since the scattering length of the interaction is of the order 

1/A^, the two-particle density jj^^ develops a short scale correlation structure on the length scale 

1/A^, which is exactly the same length scale on which the potential Vn varies. For this reason the 

two limits in (I2.16P cannot be taken independently. In order to obtain the correct Gross-Pitaevskii 

equation (jl.4p we need to take into account the correlations among the particles, and the short scale 

(2) 

structure they create in the marginal density 7)^ 

To describe the correlations among the particles we make use of the solution /Ar(x) to the zero 
energy scattering equation (j2.4p . Assuming that the function /Ar(xj — Xj) gives a good approximation 
for the correlations between particles i and j, we may expect that the one- and two-particle densities 
associated with the evolution of a condensate are given, for large but finite A^, by 



7yj(xi;x'i) ~ ipt{xi)ipt{x[) 

7^,l(a;i,X2;Xi,X2) ^ fN{xi - X2)fN{x'i - X2)(pt{xi)ipt{x2)Tpt{x[)7pt{x'2) 

Inserting this ansatz into ()2.1ip . we obtain a new self-consistent equation 



(2.17) 



-Av9t+ lim N^^ / dxf{Nx)V{Nx) \<ft\ <ft 

, N—*oo 



(2.18) 



-Acpt + 87rao\ipt\ (ft 



because of (|2.3p . This is exactly the Gross-Pitaevskii equation (II. 4p . with the correct coupling 
constant in front of the nonlinearity. 
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Note that the presence of the correlation functions /Ar(xi —X2) and fN{x'i —X2) in (j2.17p does not 
contradict complete condensation of the system at time t. On the contrary, in the weak limit — > 00, 
the function /tv converges to one, and therefore 7^"*^ and converge to \ipt){ipt\ and \ipt) {ipt\'^'^ , 
respectively. The correlations described by the function Jn can only produce nontrivial effects on 
the macroscopic dynamics of the system because of the singularity of the interaction potential V/v- 

From this heuristic argument it is clear that, in order to obtain a rigorous derivation of the Gross- 
Pitaevskii equation (j2.18p , we need to identify the short scale structure of the marginal densities and 

prove that, in a very good approximation, it can be described by the function /jv as in (j2.17p . In 

(2) 

other words, we need to show a very strong separation of scales in the marginal density 7]y ^ (and, 
more generally, in the A;-particle density 7^ j) associated with the solution of the A^-body Schrodinger 

(k) 

equation; the Gross-Pitaevskii theory can only be correct if t has a regular part, which factorizes 
for large N into the product of k copies of the orthogonal projection \(pt){(pt\, and a time independent 
singular part, due to the correlations among the particles, and described by products of the functions 

fN{xi - Xj), 1 <i,j <k. 



3 Main Results 

To prove our main results we need to assume the interaction potential to be sufficiently weak. To 
measure the strength of the potential, we introduce the dimensionless quantity 



a 



sup |x|V(x) + / p:V{x). (3.1) 



Apart from the smallness assumption on the potential, we also need to assume that the correlations 
characterizing the initial A^-particle wave function are sufficiently weak. We define therefore the 
notion of asymptotically factorized wave functions. We say that a family of permutation symmetric 
wave functions tpN is asymptotically factorized if there exists ip G L^(M'^) and, for any fixed A; > 1, 
there exists a family C/^"'^'' ^ Lg(M^(^~'^)) such that 



as N ^00. (3.2) 



It is simple to check that, if V'Af is asymptotically factorized, then it exhibits complete Bose-Einstein 
condensation in the one-particle state (p (in the sense that the one-particle density associated with 
ipN satisfy 7^'' \'p){^\ as A^ — > 00). Asymptotic factorization is therefore a stronger condition 
than complete condensation, and it provides more control on the correlations oi tpisf. 

Theorem 3.1. Assume that V{x) is a positive, smooth, spherical symmetric, and compactly sup- 
ported potential such that a (defined in f fg. 1]) ) is sufficiently small. Consider an asymptotically 
factorized family of wave functions ipN ^ L^^R^^), exhibiting complete Bose-Einstein condensation 
in a one-particle state 99 G //-"^(M^), in the sense that 

iN^ ^\^){^\ as 00 (3.3) 

where 7^^ denotes the one-particle density associated with -^at. Then, for any fixed t € the one- 
particle density 7^^^ associated with the solution il^N,t of the N -particle Schrodinger equation 112. 6\) 
satisfies 

\^t){^t\ as N ^ 00 (3.4) 
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where (ft is the solution to the time- dependent Gross- Pitaevskii equation 

idtipt = -A(/3t + 8TTao\ipt\'^(pt (3.5) 

with initial data ipt=Q = 92. 

The convergence in ()3.3p and (|3.4p is in the trace norm topology (which in this case is equivalent 
to the weak* topology defined on the space of trace class operators on M^). Moreover, from ()3.4p we 
also get convergence of higher marginal. For every /c > 1, we have 

7S-l9^t)('^tr' asiV^oo. 

Theorem 13.11 can be used to describe the dynamics of condensates satisfying the condition of 
asymptotic factorization. The following two corollaries provide examples of such initial data. 

The simplest example of A^-particle wave function satisfying the assumption of asymptotic fac- 
torization is given by a product state. 

Corollary 3.2. Under the assumptions on V{x) stated in Theorem \3.1l let V'Af(x) = Y[f=i vi^j) f^"^ 
an arbitrary 99 G i7-'^(M'^). Then, for any t G R, 

^N]t^\Vt){^t\ as N ^ 00 

where ipt is a solution of the Gross-Pitaevskii equation 13. 5\) with initial data ipt=Q = f- 

The second application of Theorem 13.11 gives a mathematical description of the results of the 
experiments depicted in the introduction. 
Let 

N N 
h'n" = Yl (-^J + ^ext(a;,)) + VN{xi - Xj) (3.6) 

j=l i<j 

with a confining potential V^xt- Let ipjsr be the ground state of H^^^^ . By [9], iI^n exhibits complete 
Bose Einstein condensation into the minimizer (/>gp of the Gross-Pitaevskii energy functional 8q-p 
defined in (jl.2p . In other words 

7^^ ^ |0GP)(0GP| asiV^cxD. 

In [5], we demonstrate that tp^ also satisfies the condition (j3.2p of asymptotic factorization. From 
this observation, we obtain the following corollary. 

Corollary 3.3. Under the assumptions on V{x) stated in Theorem \3.1{ let Tp^ he the ground state 

of ^3. and denote by 7^^^ the one-particle density associated with the solution ipN,t = e~^^'^'''tpN 
of the Schrodinger equation \2. 6]) . Then, for any fixed t G M, 

lN]t-'\'Pt){Vt\ as CO 

where (ft is the solution of the Gross-Pitaevskii equation 13. 5\) with initial data ipt=o = (pGP- 

Although the second corollary describes physically more realistic situations, also the first corollary 
has interesting consequences. In Section [21 we observed that the emergence of the scattering length in 
the Gross-Pitaevskii equation is an effect due to the correlations. The fact that the Gross-Pitaevskii 
equation describes the dynamics of the condensate also if the initial wave function is completely 



7 



uncorrelated, as in Corollary 13.21 implies that the A'"-body Schrodinger dynamics generates the 
singular correlation structure in very short times. Of course, when the wave function develops 
correlations on the length scale the energy associated with this length scale decreases; since 

the total energy is conserved by the Schrodinger evolution, we must conclude that together with 
the short scale structure at scales of order 1/N ^ the A^-body dynamics also produces oscillations on 
intermediate length scales <C ^ <C 1, which carry the excess energy (the difference between the 
energy of the factorized wave function and the energy of the wave function with correlations on the 
length scale and which have no effect on the macroscopic dynamics (because only variations 

of the wave function on length scales of order one and order affect the macroscopic dynamics 
described by the Gross-Pitaevskii equation). 



4 General Strategy of the Proof and Previous Results 

In this section we illustrate the strategy used to prove Theorem 13.11 The proof is divided into three 
main steps. 

Step 1. Compactness 0/7]^^. Recall, from ()2.7p . the definition of the marginal densities ^ni 
associated with the solution il^N,t = e^p{—iH]\ft)7p]\f of the A^-body Schrodinger equation. By defini- 
tion, for any € N and t G M, 7^^^^ is a positive operator in = JC^{L^{R^^)) (the space of trace 
class operators on L'^{M.^^)) with trace equal to one. For fixed t G M and A; > 1, it follows by standard 



general argument (Banach-Alaouglu Theorem) that the sequence {7^j}Ar>A: is compact with respect 
to the weak* topology of C\. Note here that Cj. has a weak* topology because Cj. = /C|,, where 
/Cfc = /C(L^(R^'^)) is the space of compact operators on L'^(M.^^). To make sure that we can find 

(k) 

subsequences of 7]^ ^ which converge for all times in a certain interval, we fix T > and consider the 
space C{[0,T], Cl.) of all functions of t € [0,T] with values in £]. which are continuous with respect 
to the weak* topology on Since )Ck is separable, it follows that the weak* topology on the unit 
ball of CI is metrizable; this allows us to prove the equicontinuity of the densities Xn v *o obtain 
compactness of the sequences {7^|}7V>fc in ^'([O, T], £^). 

Step 2. Convergence to an infinite hierarchy. By Step 1 we know that, as A^ —i- 00, the 
family of marginal densities Ttv,* = {7^t}fcLi least one limit point Too,* = {7^|t}fc>i in 

0j;.>^ C([0, T], >C^) with respect to the product topology. Next, we derive evolution equations for 

the limiting densities 7^\. Starting from the BBGKY hierarchy (j2.10p for the family r7v,t) we prove 
that any limit point Too,* satisfies the infinite hierarchy of equations 

k k 

idat]t = E [-^vlt]] +87rao J]Trfc+i [5{x, - Xk+i),ltf^] (4-1) 
i=i i=i 

for k>l. It is at this point, in the derivation of this infinite hierarchy, that we need to identify the 
singular part of the densities 7];^ ^ . The emergence of the scattering length in the second term on 
the right hand side of (j4.ip is due to short scale structure of 7^^^'*- 

It is worth noticing that the infinite hierarchy ()4.ip has a factorized solution. In fact, it is simple 
to see that the infinite family 

7^ = |V9i)(v.i|®'= forfc>l (4.2) 
solves (|4.ip if and only if (pt is a solution to the Gross-Pitaevskii equation (|3.5p . 
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Step 3. Uniqueness of the solution to the infinite hierarchy. To conclude the proof of Theorem [STTl 
we show that the infinite hierarchy (j4.ip has a unique solution. This implies immediately that the 

Ik) 

densities 7]y ^ converge; in fact, a compact sequence with at most one limit point is always convergent. 
Moreover, since we know that the factorized densities (j4.2p are a solution, it also follows that, for 
any k>l, 

with respect to the weak* topology of C\. 

Similar strategies have been used to obtain rigorous derivations of the nonlinear Hartree equation 

idt(pt = -^^^1 + {v *\<ft?')^t (4.3) 

for the dynamics of initially factorized wave functions in bosonic many particle mean field models, 
characterized by the Hamiltonian 

N N 

H^^ = ^-A, + -^v{x.-x,). (4.4) 

j=l i<j 

In this context, the approach outlined above was introduced by Spohn in [11], who applied it to 
derive ()4.3p in the case of a bounded potential v. In [8], Erdos and Yau extended Spohn's result to 
the case of a Coulomb interaction v{x) = ±l/|x| (partial results for the Coulomb case, in particular 
the convergence to the infinite hierarchy, were also obtained by Bardos, Golse, and Mauser, see [3|). 
More recently, Adami, Golse, and Teta used the same approach in [T] for one-dimensional systems 
with dynamics generated by a Hamiltonian of the form (j4.4p with an A^-dependent pair potential 
vn{x) = N^V{N^x), (3 < 1. In the limit — > oo, they obtain the nonlinear Schrodinger equation 

idtft = -A(/3t + bo\ft\'^(pt with bo = j V{x)dx . 

Notice that the Hamiltonian (j2.5p has the same form as the mean field Hamiltonian (j4.4p . with 
an A^-dependent pair potential vj\[{x) = N^V{Nx). Of course, one may also ask what happens if 
we consider the mean field Hamiltonian (j4.4p with the A^-dependent potential vj\[{x) = N^^V{N^x), 
for P ^ 1. If /3 < 1, the short scale structure developed by the solution of the Schrodinger equation 
is still characterized by the length scale 1/A^ (because the scattering length of N^'^^^V{Nf^x) is still 
of order 1/N); but this time the potential varies on much larger scales, of the order ^ A^^^. 
For this reason, if /? < 1, the scattering length does not appear in the effective macroscopic equation 
(Sttoo is replaced by 6o = / dx V{x)). In \Q\ (and previously in [5j for < /? < 1/2) we prove in fact 
that Corollarv 13.21 can be extended to include the case < /3 < 1 as follows. 

Theorem 4.1. Suppose iPn{'^) = Y[f=i^i^j)' f'^^ some ip G H'^(U.^). Let 'ipN,t = e~*-'^'3'^*^7v with 
the mean-field Hamiltonian 

N ^ N 

Hp,^ = E + ]v ^ Ar3/^F(Ar'3(:r, - x,)) 

j=l i<j 

for a positive, spherical symmetric, compactly supported, and smooth potential V such that a ( defined 
in 113. is sufficiently small. Let be the one-particle density associated with ipN,t- Then, if 

< /? < 1 we have, for any fixed t € M, 7^^^ \Lpt){<pt\ as N ^ 00. Here (ft is the solution to the 
nonlinear Schrodinger equation 

idtipt = -Aipt + cF\Lpt\^ipt 
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with initial data (pt=o = f o.^d with 



a 



Svroo if f3 = 1 
bo i/0</3<l 



5 Convergence to the Infinite Hierarciiy 

In this section we give some more details concerning Step 2 in the strategy outhned above. We 
consider a hmit point r^o,* = {l^t}k>i of the sequence r^v,* = {7^1}^i and we prove that T^^t 
satisfies the infinite hierarchy (j4.ip . To this end we use that, for finite A^, the family T^^t satisfies the 
BBGKY hierarchy (j2.10p . and we show the convergence of each term in (j2.10p to the corresponding 
term in the infinite hierarchy (j4.ip (the second term on the r.h.s. of (|2.10p is of smaller order and 
can be proven to vanish in the limit N — > oo). 

The main difficulty consists in proving the convergence of the last term on the right hand side 
of (I2.10p to the last term on the right hand side of (j4.ip . In particular, we need to show that in the 
limit ^ oo we can replace the potential {N — k)N'^V{N{xj — x^+i)) — N'^V{Nx) in the last term 
on the r.h.s. of (j2.10p by SiraoS^Xj — x^+i) ■ In terms of kernels we have to prove that 

dxfc+i {N^V{N{xj -Xk+i)) -8Trao6{xj - Xk+i)) 7^^^^^(xfc, x^+i, x^, x^+i) (5.1) 

as ^ oo. It is enough to prove the convergence ()5.ip in a weak sense, after testing the expression 
against a smooth fc-particle kernel j('^)(xfc; x^). Note, however, that the observable J^*^^ does not 
help to perform the integration over the variable x^+i- 

The problem here is that, formally, the A^-dependent potential N'^V{N{xj — x^+i)) does not 
converge towards 8TTaQ6{xj — Xk+i) as N ^ oo (it converges towards bo6{xj — Xk+i), with bo = 
J dxV{x)). Eq. (j5.ip is only correct because of the correlations between Xj and x^+i hidden in the 

density 7^^^^- Therefore, to prove (|5.ip . we start by factoring out the correlations explicitly, and 
by proving that, as A^ — > oo, 

^ fAjS^^fATf ^^f f AC xf 7N,t ^(xfc,Xfc+i,Xfe,Xfc+i) 

dxk+i [N V{N{xj - Xk+i)jfN{xj - Xk+i) - 8TTaod{xj - Xk+i)) r > 0, 

jN\Xj — a^fc+ij 

(5.2) 

where fN{x) is the solution to the zero energy scattering equation (j2.4p . Then, in a second step, we use 
the fact that /at — > 1 in the weak limit A^ — > oo, to prove that the ratio 7]y ^ / fis[{xj —Xk+i) converges 

to the same limiting density '^oot as 7]^^ . Eq. (|5.2I) looks now much better than (|5.ip because, 
formally, N'^V{N{xj — Xk+i))fN{xj — x^+i) does converge to 87rao(5(xj — x^+i). To prove that (|5.2p 

is indeed correct, we only need some regularity of the ratio 7^^^''(xfc, x^+i; x'^, XA:+i)//Ar(xj — xt+i) 
in the variables Xj and Xk+i- In terms of the A^-particle wave function tpN,t we need regularity of 
ipN,t{^)/ fN{xi — Xj) in the variables Xi, Xj, for any i ^ j- To establish the required regularity we use 
the following energy estimate. 

Proposition 5.1. Consider the Hamiltonian H]\j defined in i2.5\) . with a positive, spherical sym- 
metric, smooth and compactly supported potential V . Suppose that a (defined in \3. is sufficiently 
small. Then there exists C = C{a) > such that 



{iP, Hl^) > CN^ / dx 
for all i / j and for all G Ll{R^^ , dx). 



fN{Xi - 



(5.3) 
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Making use of this energy estimate it is possible to deduce strong a-priori bounds on the solution 
V'Art of the Schrodinger equation (j2.6p . These bounds have the form 



/ 



dx 



jNyXi — Xj 



< C (5.4) 



uniformly in iV G N and i G R. To prove (j5.4p we use that, by (j5.3p . and because of the conservation 
of the energy along the time evolution, 



dx 



jN{Xi - Xj) 



2 



< CN~-'{il^M,uHl^l,N,t) = CN-^{%bN,o,Hi,iPN,o) ■ (5.5) 



From ()5.5p and using an approximation argument on the initial wave function to make sure that the 
expectation of Hfj^ at time t = is of the order N'^, we obtain (j5.4p . 

The bounds (j5.4p are then sufficient to prove the convergence (|5.ip (using a non-standard Poincare 
inequality; see Lemma 7.2 in [6]). 

Remark that the a-priori bounds (j5.4p do not hold true if we do not divide the solution ipN,t of 
the Schrodinger equation by fN{xi — xj) (replacing '4'N,ti^)/ fwixi — Xj) by V'Ar(x) the integral in 
(j5.4p would be of order A''). It is only after removing the singular factor fN{xi — Xj) from iljN,t{?^) 
that we can prove useful bounds on the regular part of the wave function. 

It is through the a-priori bounds (j5.4p that we identify the correlation structure of the wave 
function 'ij^N^t and that we show that, when Xi and Xj are close to each other, ipi\i^t{^) can be 
approximated by the time independent singular factor fiy{xi — Xj), which varies on the length scale 
1/N, multiplied with a regular part (regular in the sense that it satisfy the bounds (|5.4p ). It is 
therefore through (15. 4[) that we establish the strong separation of scales in the wave function ^l^N,t 

(k) 

and in the marginal densities 7]^ ^ which is of fundamental importance for the Gross-Pitaevskii theory. 

Since it is quite short and it shows why the solution fNixi — xj) to the zero energy scattering 
equation ()2.ip can be used to describe the two-particle correlations, we reproduce in the following 
the proof Proposition l5.ll Note that this is the only step in the proof of our main theorem where the 
smallness of constant a, measuring the strength of the interaction potential, is used. The positivity 
of the interaction potential, on the other hand, also plays an important role in many other parts of 
the proof. 

Proof of Proposition 15. ii We decompose the Hamiltonian (|2.5p as 

^ 1 
Hn = ^ hj with hj = —Aj + - ^ VNixi — Xj) . 

For an arbitrary permutation symmetric wave function ip and for any fixed i 7^ j, we have 

{iP,H^NiP) = N{i;,h^iP) + N{N - l){iP,hihjiP) > NiN - 1)(V^, /ii/i.-^) . 
Using the positivity of the potential, we find 

(V, H%i;) > N{N - 1) (i^, (^-A, + ^Nixi - xj)^ (^-Aj + Wn{x^ - x,)^ . (5.6) 

Next, we define by ip{x) = /Ar(xj — Xj) (p{x) {(f) is well defined because fN{x) > for all x £ 
note that the definition of the function (/> depends on the choice of i,j. Then 

1 A ™ _ A , (.'^fN){xi-Xj)^^^^^ ^ V/Ar(xi-Xj), 



Ai ifNixi - x,)0(x)) = A,,/.(x) + ^ /"/^ ^(/.(x) + f Vi</)(x) 



/^(X,-X,-) ' ^■'^'^ ' ''^' '^ fN{Xi-X,) -^^"^ ' fN{x. 
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From (pJ]) it follows that 

1 



fN{Xi - Xj) 

and analogously 

1 



fN{Xi - Xj) 

where we defined 



-Ai + ^VNiXi - Xj)^ fN{Xi - Xj)4>{-K) = Li(j){x) 
-Aj + ^VN{Xi - Xj)^ fN{Xi - Xj)(t>{x) = Lj4>{x) 



Li, = -Ai + 2—— r^V^ for £ = i,j. 

jN{Xi - Xj) 

Remark that, for i = the operator satisfies 

j dx/^(xi-Xj) Li(j){x) '0(x) = j dxfff{xi-Xj) (j){x) Li^{y.) = J dx /^(xj-Xj) VeH^) V^V^W • 



Therefore, from (j5.6p . we obtain 

(V-, H^N^) > N{N - 1) / dx f^{xi - Xj) Li^(x) Lj <^(x) 



= N{N - 1) / dx fUxi - Xj) VicPi^) ViLj (/>(x) 

= N{N - 1) /" dx /^(xi - Xj) Vi0(x) Lj V,(/>(x) 

^ . _ (5.7) 

+ N{N - 1) y dx - Xj) Vi(/.(x) [V„ Lj](/.(x) 

= Af(Af-l) y dx/2,(xi-x,-) |V,Vi</)(x)|2 

+ iV(iV-l) / dx/2,(x,-x,) ("v.^M^i^) V,0(x)V,</<(x). 

J V jN\Xi — Xj) J 

To control the second term on the right hand side of the last equation we use bounds on the function 
/at, which can be derived from the zero energy scattering equation (j2.ip : 



1 - Ca < /iv(x) < 1, lV/jv(x)| <C^, |v2/iv(x)| (5.8) 

|x| |x|^ 

for constants C independent of N and of the potential V (recall the definition of the dimensionless 
constant a from (j3.ip ). Therefore, for a < 1, 

/dx/^(x.-x,) fv.^M^i^) V.0(x) V,,^(x) 

J V /Af(a^i - Xj) J 

<Ca [dx ^—-^ |V,(A(x)| \Vj^{x)\ 

J 1*^2 •^J I 

<Ca [dx ^ (|Vi</.(x)|2 + lV,(A(x)|2) 

J I X j X J I 

<Cajdx |V,Vj0(x)|2 

where we used Hardy inequality. Thus, from (j5.7p . and using again the first bound in (j5.8p . we obtain 

{iP,H%iP)> N{N-l){l-Ca) J dx\ViVj(P{x)f 
which implies (15. 3p . □ 



(5.9) 
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6 Uniqueness of the Solution to the Infinite Hierarchy 

In this section we discuss the main ideas used to prove the uniqueness of the solution to the infinite 
hierarchy (Step 3 in the strategy outhned in Section H]). 

First of all, we need to specify in which class of family of densities Fj = {7^ '^^}a;>i we want to prove 
the uniqueness of the solution to the infinite hierarchy (j4.ip . Clearly, the proof of the uniqueness 
is simpler if we can restrict our attention to smaller classes. But of course, in order to apply the 
uniqueness result to prove Theorem 13.11 we need to make sure that any limit point of the sequence 
TN,t = {T^tlfcLi is ™ the class for which we can prove uniqueness. 

We are going to prove uniqueness for all families Fj = {7t'^^}fc>i € C([0, T], >C^) with 



17^11... :=Tr 



(1 - A.y/^ ... (1 - A,)V2 _ A,) V2 . . . (1 _ A.y/^ < (6.1 



for all t € [0, T] and for all /c > 1 (with a constant C independent of k). 

The following proposition guarantees that any limit point of the sequence F^v,* satisfies ()6.1 



Proposition 6.1. Assume the same conditions as in Proposition \5.1l Suppose that Foo,t = {7,^\}fc>i 
is a limit point ofT^^t = {7/^j}fcLi with respect to the product topology on 0^>;^ C{[{),T], C\). Then 

(k) ~ 

loo t — ^ there exists a constant C such that 

Tr(l-Ai)...(l-Afc)7S<C'= (6.2) 

for allk>l and t G [0,T]. 

Because of Proposition 16.11 it is enough to prove the uniqueness of the infinite hierarchy (j4.ip in 
the following sense. 

Theorem 6.2. Suppose that F = {7'''^'*}fc>i is such that 

h^'^llw. <C'= (6.3) 

for all k > 1 (the norm W-Wuk ^-^ defined in i6.1\) ). Then there exists at most one solution F^ = 
{7f ^}fc>i G ©C([0,r], A) of O such that Ft=o = T and 

\\iP\\h,<C' (6.4) 

for all k > 1 and all t G [0, T] (with the same constant C as in \6. 

In the next two subsections we explain the main ideas of the proofs of Proposition 16.11 and 
Theorem [ 

6.1 Higher Order Energy Estimates 

The main difficulty in proving Proposition 16.11 is the fact that the estimate (16. 2p does not hold true 

(k) (k) 

if we replace 7^0 j by the marginal density Tn t- More precisely, 

Tr (1 - Ai) ... (1 - Akh^j^l ^ ^'^ (6.5) 
cannot hold true with a constant C independent of N. In fact, for finite and k > 1, the /c-particle 

(k) 

density jj^^ still contains the short scale structure due to the correlations among the particles. 
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Therefore, when we take derivatives of as in (j6.5p . the singular structure (which varies on a 
length scale of order generates contributions which diverge in the limit — > oo. 

To overcome this problem, we cutoff the wave function ipN,t when two or more particles come 
at distances smaller than some intermediate length scale £, with A^~^ <C ^ «C 1 (more precisely, the 
cutoff will be effective only when one or more particles come close to one of the variable Xj over 
which we want to take derivatives). For fixed j = 1, . . . , A^, we define 6j G C°°(M^^) such that 

Q ( \ i ^ ~ Xj\ ^ i for all i j 

■' \ if there exists i ^ j with |xj — Xj\ < i 

It is important, for our analysis, that 9j controls its derivatives (in the sense that, for example, 
|Vj0j| < Cl^^Oy^); for this reason we cannot use standard compactly supported cutoffs, but instead 
we have to construct appropriate functions which decay exponentially when particles come close 
together. Making use of the functions Oj{x.), we prove the following higher order energy estimates. 

Proposition 6.3. Choose £ <^ 1 such that Nl'^ ^ 1. Suppose that a is small enough. Then there 
exist constants Ci and C2 such that, for any ip € L^(M^^), 

{iP, {Hn + CiN)'iP) > C2N'' J dx 01 (x) . . . 0fc_i(x) |Vi . . . VfcV(x)P . (6.6) 

The meaning of the bounds (16. 6p is clear. We can control the L^-norm of the k-th derivative 
Vi . . . VfcV by the expectation of the k-th. power of the energy per particle, if we only integrate over 
configurations where the first k — 1 particles are "isolated" (in the sense that there is no particle at 
distances smaller than i from xi,X2, ■ ■ ■ ,Xk-i)- In this sense the energy estimate in Proposition 15.11 
(which, compared with Proposition 16.31 is restricted to /c = 2) is more precise than ()6.6p . because 
it identifies and controls the singularity of the wave function exactly in the region cutoff from the 
integral on the right side of (j6.6p . The point is that, while Proposition 15.11 is used to identify the 

(k) 

two-particle correlations in the marginal densities 'jj^^ (which are essential for the emergence of 
the scattering length oq in the infinite hierarchy (j4.ip l. we only need Proposition 16.31 to establish 
properties of the limiting densities; this is why we can introduce cutoffs in ()6.6p . provided we can 
show their effect to vanish in the limit N ^ 00. 

Note that we can allow one "free derivative"; in (j6.6p we take the derivative over although 
there is no cutoff 0a:(x). The reason is that the correlation structure becomes singular, in the L^- 
sense, only when we derive it twice (if one uses the zero energy solution /tv introduced in (j2.ip to 
describe the correlations, this can be seen by observing that V/Ar(2;) ~ which is locally square 

integrable). Remark that the condition Ni'^ ^ 1 is a consequence of the fact that, if i is too small, 
the error due to the localization of the kinetic energy on distances of order i cannot be controlled. 
The proof of Proposition 16.31 is based on induction over k; for details see Section 9 in [6]. 

From the estimates ()6.6p . using the preservation of the expectation of along the time evolution 
and a regularization of the initial A^-particle wave function ^pN, we obtain the following bounds for 
the solution ipN,t = e~^^'^^'ipN of the Schrodinger equation (j2.6p . 

J dx 01 (x) . . . 0fc_i(x) [Vi . . . VfcVjv,t(x)|' < (6.7) 

uniformly in and t, and for all /c > 1. Translating these bounds in the language of the density 
matrix 'yN,t, we obtain 

TV0i...0fc_iVi...Vfc77v,tVl...V^ (6.8) 
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The idea now is to use the freedom in the choice of the cutoff length L If we fix the position of all 
particles but Xj, it is clear that the cutoff 9j is effective in a volume at most of the order Ni"^. If 
we choose now I such that Nl'^ — > as ^ oo (which is of course compatible with the condition 
that » 1), then we can expect that, in the limit of large A^, the cutoff becomes negligible. This 
approach yields in fact the desired results; starting from (|6.8|) . and choosing I such that A^^^ <C 1, 
we can complete the proof of Proposition 16.11 (see Proposition 6.3 in [6] for more details). 



6.2 Expansion in Feynman Graphs 

To prove Theorem 16. 2^ we start by rewriting the infinite hierarchy ()4.ip in the integral form 
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(6.9) 



= ZY(^')(t)7o+ / dsZYW(t-s)S('=)7f+^), 
Jo 

where lA^^\t) denotes the free evolution of k particles, 

and the collision operator B^^^ maps {k + l)-particle operators into /c-particle operators according to 



i=i 

(recall that Trfc+i denotes the partial trace over the {k + l)-th particle). 
Iterating (j6.9p n times we obtain the Duhamel type series 

n-l 

7f =Z^('=)(t)7f + E^S + ^S 



(6.10) 



(6.11) 



m=l 



with 



dSm U^'Ht - Si)i?(%('=+^)(S1 - S2)B('+'^ . . . S(^+'"-l)Z^('=+-)(s^)7('=+™) 

- Xk+l), 

U^'^'-'Hsi - S2)Trfc+2 \6{XJ, - Xk+2), ■ ■ ■ TTk+m U{X,^ - Xk+m)M''''"'Hsmhjf^"'^ ■ ■ 



k k+1 k+m „i 

jl = lj2 = l jm = l-^" 



(6.12) 



and the error term 



(fe) 



t f-Sl 

dsi / dso ■ ■ ■ 



^0 



Sn-l 



dsn U^'^Ht - si)B^^^U^^+^\si - S2)B^''+^^ . . . . 



(6.13) 

Note that the error term (j6.13p has exactly the same form as the terms in (|6.12p . with the only 
difference that the last free evolution is replaced by the full evolution 75^"''™^ . 
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Vertices: 



2k roots 




Figure 1: A Feynman graph in J^rn,k and its two types of vertices 

To prove the uniqueness of the infinite hierarchy, it is enough to prove that the error term (j6.13p 
converges to zero as n — > oo (in some norm, or even only after testing it against a sufficiently large 
class of smooth observables) . The main problem here is that the delta function in the collision 
operator B'^^^ cannot be controlled by the kinetic energy (in the sense that, in three dimensions, the 
operator inequality 5{x) < C(l — A) does not hold true). For this reason, the a-priori estimates 
llTt'^^IIWfe — sufficient to show that (j6.13p converges to zero, as n — > oo. Instead, we have 

to make use of the smoothing effects of the free evolutions lA^^~^^\sj — Sj+i) in (in a similar 

way, Stricharzt estimates are used to prove the well-posedness of nonlinear Schrodinger equations). 
To this end, we rewrite each term in the series (j6.1ip as a sum of contributions associated with certain 
Feynman graphs, and then we prove the convergence of the Duhamel expansion by controlling each 
contribution separately. 

The details of the diagrammatic expansion can be found in Section 9 of [5] . Here we only present 
the main ideas. We start by considering the term sj^^ in (j6.12p . After multiplying it with a compact 
fc-particle observable J^^'^ and taking the trace, we expand the result as 

TVJ«d'i= E ^A.* (6.14) 

where K\^t is the contribution associated with the Feynman graph A. Here Trn,k denotes the set of 
all graphs consisting of 2k disjoint, paired, oriented, and rooted trees with m vertices. An example 
of a graph in J-m,k is drawn in Figure [H Each vertex has one of the two forms drawn in Figure [U 
with one "father" -edge on the left (closer to the root of the tree) and three "son" -edges on the right. 
One of the son edge is marked (the one drawn on the same level as the father edge; the other two 
son edges are drawn below). Graphs in J-m,k have 2k + 3m edges, 2k roots (the edges on the very 
left), and 2k + 2m leaves (the edges on the very right). It is possible to show that the number of 
different graphs in Tm,k is bounded by 2^'""'"'^. 

The particular form of the graphs in J-m.k is due to the quantum mechanical nature of the 
expansion; the presence of a commutator in the collision operator ()6.10p implies that, for every 
we can choose whether to write the interaction on the left or on the right of the 
density. When we draw the corresponding vertex in a graph in J-'m,k, we have to choose whether to 
attach it on the incoming or on the outgoing edge. 
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Graphs in J-m,k are characterized by a natural partial ordering among the vertices [v -< v' if 
the vertex v is on the path from v' to the roots); there is, however, no total ordering. The absence 
of total ordering among the vertices is the consequence of a rearrangement of the summands on 
the r.h.s. of (j6.12p : by removing the order between times associated with non-ordered vertices we 
significantly reduce the number of terms in the expansion. In fact, while (|6.12|) contains (m + ky./kl 
summands, in (j6.14p we are only summing over 2^™+'^ contributions. The price we have to pay is 
that the apparent gain of a factor 1/m! because of the ordering of the time integrals in (j6.12p is lost 
in the new expansion ()6.14p . However, since the time integrations are already needed to smooth out 
singularities, and since they cannot be used simultaneously for smoothing and for gaining a factor 
1/m!, it seems very difficult to make use of the apparent factor 1/m! in (j6.12p . In fact, we find that 
the expansion (j6.14p is better suited for analyzing the cumulative space-time smoothing effects of 
the multiple free evolutions than (j6.12p . 

Because of the pairing of the 2k trees, there is a natural pairing between the 2k roots of the 
graph. Moreover, it is also possible to define a natural pairing of the leaves of the graph (this is 
evident in Figure [1]) ; two leaves ii and £2 are paired if there exists an edge ei on the path from ii 
back to the roots, and an edge 62 on the path from £2 to the roots, such that ei and 62 are the two 
unmarked son-edges of the same vertex (or, if there is no unmarked sons in the path from ii and £2 
to the roots, if the two roots connected to ii and £2 are paired). 

For A € J^m.kj we denote by E{A), V{A), i?(A) and L{A) the set of all edges, vertices, roots 
and, respectively, leaves in the graph A. For every edge e S E{A), we introduce a three-dimensional 
momentum variable Pe and a one-dimensional frequency variable Qg. Then, denoting by ^y^'^~^™'^ and 
by J^^^ the kernels of the density ^^'^^"^^ and of the observable J^^'^ in Fourier space, the contribution 
K\^t in (|6.14p is given by 

eGE(A) ^ DGV(A) Veef / \eev / 

X exp I -it ^ Te{ae+iTeHe) 1 J^''^ {{Pe}eeR(A)) ^q''^"^^ {{Pe}eeL(A)) ■ 
\ eeH(A) / 

(6.15) 

Here Tg = ±1, according to the orientation of the edge e. We observe from (j6.15p that the momenta 
of the roots of A are the variables of the kernel of J^^^ , while the momenta of the leaves of A are the 
variables of the kernel of -y^^~^"^^ (this also explain why roots and leaves of A need to be paired). 

The denominators {ae—p1+iTei-ie)~^ are called propagators; they correspond to the free evolutions 
in the expansion (j6.12p and they enter the expression (|6.15p through the formula 

da 15 

a — p'' + i/j. 

(here and in (j6.15p the measure da is defined by da = d'a/(27ri) where d'a is the Lebesgue measure 
on M). 

The regularization factors /ig in (j6.15p have to be chosen such that /Xfather = Yle= son /^e every 
vertex. The delta-functions in (j6.15p express momentum and frequency conservation (the sum over 
e & V denotes the sum over all edges adjacent to the vertex v; here ±ae = ae if the edge points 
towards the vertex, while ztag = — ag if the edge points out of the vertex, and analogously for ±Pe)- 
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An analogous expansion can be obtained for the error term -ql^] in (j6.13p . The problem now 
is to analyze the integral (|6.15p (and the corresponding integral for the error term). Through an 
appropriate choice of the regularization factors //g one can extract the time dependence of K\^t and 
show that 

\K,A<c^--t-/' j n n 

■' e&E(V) -^^^ v&v{r) \eev ) \eev j (6.16) 

X p'^) {{Ve^emV)) I ItS'^"^ ({Pe}e.L(r)) | 

where we introduced the notation (x) = (1 + x^)"*^/^. 

Because of the singularity of the interaction at zero, we may be faced here with an ultraviolet 
problem; we have to show that all integrations in (j6.16p are finite in the regime of large momenta 

and large frequency. Because of (j6.3p . we know that the kernel 7o'^~'''^^({Pe}eeL(A)) i'^ (|6.16p provides 
decay in the momenta of the leaves. From (j6.3p we have, in momentum space, 



y dpi . . . dp, (p? + 1) . . . {Vl + 1) lt\vx. ■ ■ ■ ,Pn;Pl, . . . ,Pn) < 

for all n > 1. Power counting implies that 

\lt"'\{Pe}eeLiA))\ < U {Pe)-"' ■ (6.17) 

eGL(A) 

Using this decay in the momenta of the leaves and the decay of the propagators (ag— Pg)~^, e € E{A), 
we can prove the finiteness of all the momentum and frequency integrals in (j6.15p . Heuristically, this 
can be seen using a simple power counting argument. Fix k ^ 1, and cutoff all momenta |pe| ^ and 
all frequencies \ae\ > k^. Each pe-integral scales then as k^, and each ag-integral scales as k^. Since 
we have 2k + 3m edges in A, we have 2k + 3m momentum- and frequency integrations. However, 
because of the m delta functions (due to momentum and frequency conservation), we effectively only 
have to perform 2k + 2m momentum- and frequency-integrations. Therefore the whole integral in 
()6.15p carries a volume factor of the order ^5(2fc+2m) _ ^lofc+iom^ Now, since there are 2k + 2m 
leaves in the graph A, the estimate (j6.17p guarantees a decay of the order ^-5/2(2fc+2m) _ ^-5/c-5m^ 
The 2k + 3m propagators, on the other hand, provide a decay of the order /--2(2/c+3m) _ ^-ik-Gm^ 
Choosing the observable J^'^^ so that J*^'^^ decays sufficiently fast at infinity, we can also gain an 
additional decay k~^'^. Since 

,^10fc+10m ~5fc— 5m— 4fc— 6m— 6fc ^—m—5k ^ -i 

rv ' rv — ru 1 

for /{ ^ 1, we can expect (|6.15p to converge in the large momentum and large frequency regime. 
Remark the importance of the decay provided by the free evolution (through the propagators); 
without making use of it, we would not be able to prove the uniqueness of the infinite hierarchy. 

This heuristic argument is clearly far from rigorous. To obtain a rigorous proof, we use an 
integration scheme dictated by the structure of the graph A; we start by integrating the momenta 
and the frequency of the leaves (for which (j6.17p provides sufficient decay). The point here is that 
when we perform the integrations over the momenta of the leaves we have to propagate the decay 
to the next edges on the left. We move iteratively from the right to the left of the graph, until we 
reach the roots; at every step we integrate the frequencies and momenta of the son edges of a fixed 
vertex and as a result we obtain decay in the momentum of the father edge. When we reach the 
roots, we use the decay of the kernel J^^^ to complete the integration scheme. In a typical step, we 
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P„ a. 



Pw«v 



P, a, 



Figure 2: Integration scheme: a typical vertex 

consider a vertex as the one drawn in Figure [2] and we assume to have decay in the momenta of the 
three son-edges, in the form |pe|~^; e = u,d,w (for some 2 < A < 5/2). Then we integrate over 
the frequencies au,ad,aw and the momenta Pu,Pd,Pw of the son-edges and as a result we obtain a 
decaying factor |pr|~^ in the momentum of the father edge. In other words, we prove bounds of the 
form 

daudaddawdpudpddpy; S{ar = au + ad - aw)S{pr =Pu+Pd- Vw) ^ const 
\Pu\^\pd\^\Pw\^ {au-pl){ad-pi){aw-pi) ~ \Vr\^ ' 

Power counting implies that (j6.18p can only be correct if A > 2. On the other hand, to start the 
integration scheme we need A < 5/2 (from (I6.17P this is the decay in the momenta of the leaves, 
obtained from the a-priori estimates). It turns out that, choosing A = 2 -|- e for a sufficiently 
small e > 0, (j6.18p can be made precise, and the integration scheme can be completed. 

After integrating all the frequency and momentum variables, from (I6.16|) we obtain that 

for every A G J'm.k- Since the number of diagrams in J-m,k is bounded by (7*^+™, it follows immediately 
that 

Note that, from (j6.12p . one may expect S^^^ to be proportional to t"^. The reason why we only get 
a bound proportional to t^/^ is that we efi^ectively use part of the time integration to control the 
singularity of the potentials. 

Note that the only property of used in the analysis of (j6.15p is the estimate (|6.3p , which 

provides the necessary decay in the momenta of the leaves. However, since the a-priori bound (16. 4p 
hold uniformly in time, we can use a similar argument to bound the contribution arising from the 
error term rl"^\ in (|6.13p (as explained above, also r^-^l can be expanded analogously to (|6.14p . with 
contributions associated to Feynman graphs similar to (16.15P : the difference, of course, is that these 
contributions will depend on 7!'^'''"^ for all s G [0, t], while (j6.15p only depends on the initial data). 
Thus, we also obtain 

TrjC^^T/g <C'=+"t"/^ (6.19) 

This bound immediately implies the uniqueness. In fact, given two solutions Fi^t = {'^^t}k>i and 

^2,t = {l^t}k>i of the infinite hierarchy ()4.ip . both satisfying the a-priori bounds (|6.4p and with 
the same initial data, we can expand both in a Duhamel series of order n as in ()6.1ip . If we fix 
k > 1, and consider the difference between 7];/ and 72/, all terms (|6.12p cancel out because they 
only depend on the initial data. Therefore, from (I6.19p . we immediately obtain that, for arbitrary 
(sufficiently smooth) compact /c-particle operators 
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Since it is independent of n, the left side has to vanish for all t < This proves uniqueness 

for short times. But then, since the a-priori bounds hold uniformly in time, the argument can be 
repeated to prove uniqueness for all times. 
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